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ABSTRACT 

We present the results of a long-term monitoring campaign of the Ha line in a sample 
of bright OB-supergiants (07.5-B9) that aims at detecting rotationally modulated 
changes potentially related to the existence of large-scale wind structures. A total of 
22 objects were monitored during 36 nights spread over 6 months in 2001-2002. Co- 
ordinated broad-band photometric observations were also obtained for some targets. 
Conspicuous evidence for variability in Ha is found for the stars displaying a fea- 
ture contaminated by wind emission. Most changes take place on a daily time-scale, 
although hourly variations are also occasionally detected. Convincing evidence for a 
cyclical pattern of variability in Ha has been found in 2 stars: HD 14134 and HD 42087 
(periodic signals are also detected in other stars, but independent confirmation is re- 
quired) . Rotational modulation is suggested from the similarity between the observed 
recurrence time-scales (in the range 13-25 days) and estimated periods of stellar rota- 
tion. We call attention to the atypical case of HD 14134 which exhibits a clear 12.8-d 
periodicity both in the photometric and in the spectroscopic data sets. This places 
this object among a handful of early-type stars where one may observe a clear link 
between extended wind structures and photospheric disturbances. Further modelling 
may test the hypothesis that azimuthally-extended wind streams are responsible for 
the patterns of spectral variability in our target stars. 

Key words: stars: early-type - stars: supergiants - stars: emission-line - stars: winds, 
outflows - stars: rotation. 



1 INTRODUCTION 

One of the major legacies of the International Ultraviolet Ex- 
plorer (IUE) was to establish the ubiquity of UV line-profile 
variability in OB stars. The pattern of variability predom- 
inantly consists of extra absorption features (the so-called 
'discrete absorption components'; hereafter DACs) moving 
blueward across the (unsaturated) absorption troughs of P 
Cygni lines (e.g., Kaper et al. 1996). The repetitive nature 
of DACs and the fact that the recurrence time-scale is com- 
mensurate with the (otherwise ill-defined) stellar rotational 
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period suggest that rotation plays a key role in inducing the 
patterns of variability (e.g., Kaper et al. 1999). UV time- 
series with a long temporal baseline (i.e., sampling several 
rotational periods) have revealed that OB stars exhibit a 
great diversity of behaviours that cannot always be strictly 
related to the stellar rotation period (e.g., Howarth, Prinja & 
Massa 1995; Prinja, Massa & Fullerton 1995; Prinja, Massa 
& Fullerton 2002) . Nevertheless, this global phenomenon has 
been successfully interpreted in several cases in the frame- 
work of a model in which the changes are caused by the 
rotational modulation of large-scale, azimuthally-extended 
wind streams (e.g., Fullerton et al. 1997). 

The hydrodynamical simulations of Cranmer & Owocki 
(1996) suggest that large-scale photospheric perturbations 
(e.g., bright 'spots' causing a local enhancement of the ra- 
diative driving force) may lead to the development of these 
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structures. Although the underlying physical mechanism 
remains elusive, pulsational instability can be invoked as 
a trigger of the wind variability (e.g., Fullerton, Gies & 
Bolton 1996). Large-scale or 'surface' magnetic field struc- 
tures could be a viable alternative in view of the growing evi- 
dence for (weak) magnetic fields in OB stars (e.g., Henrichs, 
Neiner & Geers 2003). Photospheric perturbations arising 
from these two phenomena are likely to have a significant 
impact on the global wind morphology (e.g., ud-Doula & 
Owocki 2002). 

With the demise of IUE, further progress on this issue 
is likely to stem in the foreseeable future from time-resolved 
studies of line variability in the optical. Firstly, intensive UV 
observations have been carried out over several rotation pe- 
riods for few stars (Massa et al. f995; Prinja et al. 1998; de 
Jong et al. 2001). For practical reasons, an optical survey is 
much more amenable at collecting high-quality observations 
covering several cycles for a large, representative sample of 
OB stars. Secondly, in contrast to the UV resonance lines, 
optical recombination lines are (at least partly) formed in 
OB supergiants in the dense, strongly accelerating part of 
the outflow close to the stellar core. In virtue of their density- 
squared dependence, they are good probes of the large den- 
sity gradients expected to prevail in the co-rotating interac- 
tion regions (Cranmer & Owocki 1996). Revealing a cyclical 
pattern of variability in these lines would give further cre- 
dence to the idea that the origin of anisotropic outflows in 
OB stars is deep-seated and causally linked to processes tak- 
ing place at the stellar photosphere. Here we discuss the tem- 
poral behaviour of Ha, which is typically formed at r~ 1.5 
i?* in O-type stars (e.g., Prinja, Fullerton & Crowther 1996). 
Although the line-formation region is likely to be more ex- 
tended in B-type supergiants, this transition still probes the 
few inner radii of the stellar wind. 

Previous surveys have been very useful in documenting 
Ha line-profile variability in OB stars (e.g., Ebbets 1982), 
but frequently suffered from a poor temporal sampling ham- 
pering the study of the line-profile variations on a rotational 
time-scale (up to several weeks in the case of B supergiants). 
Despite the persistent and numerous searches for optical 
line-profile variability in OB stars, evidence for rotational 
modulation has been found in only few cases (e.g., Moffat 
& Michaud 1981; Kaper et al. 1997, 1998). In the UV and 
optical domains the time-scales associated with the changes 
are often identical, within the uncertainties (e.g., Howarth 
et al. 1995; Kaper et al. 1999; de Jong et al. 2001; Kaufer, 
Prinja & Stahl 2002). This suggests that Ha variability and 
DACs diagnose the same underlying physical phenomenon. 

The long-term monitoring of the Ha line presented in 
this paper aims at better establishing the incidence of rota- 
tionaly modulated variability in OB supergiants. To poten- 
tially link the photospheric and wind activity, coordinated 
broad-band photometic observations were also carried out 
for a subset of this sample. Our survey primarily targets 
early B-type stars and can thus be regarded as comple- 
mentary to others that concentrate on O-type (Kaper 1998; 
Kaper et al. 1998) or late B- to early A-type supergiants 
(Kaufer et al. 1996). The results presented in the following 
supersede previous preliminary reports (Morel et al. 2004a, 
b). 



2 TARGET SELECTION AND TEMPORAL 
SAMPLING OF THE OBSERVATIONS 

Our targets are drawn from a magnitude-limited (V<7.5) 
sample of early-type supergiants (07.5-B9) with evidence 
for an emission-like Ha profile or, in the case of stars with 
an absorption feature, for previous claims of a wind-related 
variability (e.g., filling-in of the profile). Known close bina- 
ries were discarded from this list in order to avoid contam- 
ination from the wind collision effects (e.g., Thaller et al. 
2001). The wide binary orbit of HD 37742 rules out the for- 
mation of a bow shock created from the interaction with an 
O-type companion (Hummel et al. 2000). In our data only 
HD 13854 and HD 47240 present significant velocity varia- 
tions of the photospheric He I A6678.15 line which might be 
attributed to binary motion (see below). Dwarfs and giant 
stars were not included in this list, as Ha is of photospheric 
origin in the vast majority of these objects. The only ex- 
ception was the 07.5 giant HD 24912 (£ Per) because of 
previous claims of a cyclical, wind-related pattern of vari- 
ability in Ha (Kaper et al. 1997; de Jong et al. 2001). For 
the other stars in our list, no clear evidence for cyclical line- 
profile variations in optical wind lines was found prior to our 
survey. Some physical parameters of our program stars are 
given in Table 1. 

The basic requirement of this survey was to achieve a 
time sampling allowing a proper assessment of the variabil- 
ity pattern on a rotational time-scale (typically 1-2 weeks). 
Although the rotational periods of our program stars are un- 
known, the critical and projected rotational velocities can be 
used, along with the estimates of the stellar radii, to set a 
lower and an upper limit on this quantity. As can be seen in 
Table 2, the estimated rotational periods span a wide range 
of values (1.2-144 d). Our 22 program stars were monitored 
during a total of 36 nights spread over 6 months, with 1-2 
spectra obtained every night. This long-term monitoring was 
necessary to cover at least one rotational cycle for the slowly 
rotating objects. We also obtained intranight observations, 
expecting significant hourly changes in some cases (e.g., HD 
24912). In general, this temporal sampling appears to be 
adequate for our purpose (see Table 3). For 5 stars, how- 
ever, the observations were too sparse to reliably investigate 
a potential cyclical behaviour (HD 14818, HD 21291, HD 
119608, HD 151804, and HD 152236). No period search was 
performed in this case (see Section 5.2), but a qualitative 
description of the changes is presented in Section 4. 



3 OBSERVATIONS AND REDUCTION 
PROCEDURE 

3.1 Spectroscopy 

Our observations were obtained during 4 observing runs con- 
ducted between 2001 November and 2002 May at the 40-inch 
telescope of the Vainu Bappu Observatory (Kavalur, India). 
The spectograph was equipped with a TK lk CCD chip. 
We used the 651 1/mm grating blazed in the 1st order at 
7500 A, which yields at Ha a reciprocal dispersion of 1.45 
A pixel -1 (we used a 120 /im-wide slit), thus covering the 
spectral region 5810-7205 A. The slit of the spectrograph 
was rotated for HD 37742 in order to avoid contamination 
from the visual companion. The speckle interferometric ob- 
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Table 1. Physical parameters of the program stars: spectral types and wind terminal velocities, Voo, from Howarth et al. (1997) and 
references therein; effective temperatures from the calibration of Humphreys & MacElroy (1984); stellar luminosities calculated from the 
calibrated My values and bolomctric corrections of Humphreys & MacElroy (1984); stellar masses from Lamers (1981) (the numbers in 
brackets are estimates based on his values for similar spectral types and luminosity classes). This source was used because it constitutes 
the most homogeneous database, but we warn the reader that these values may be significantly revised (see, e.g., Herrero, Puis & Najarro 
2002); stellar radii determined from the effective temperatures and luminosities; mass-loss rates, M, calculated from the empirical relation 
for galactic OBA supergiants of Lamers & Cassinclli (1996) (the value for HD 24912 was derived from their relation for giants); projected 
rotational velocities, tisini, from Howarth et al. (1997); critical rotational velocities, v CI n, calculated following Howarth & Prinja (1989) 
(the values in brackets arc poorly determined because of an uncertain stellar mass). 



Name 


V 


Spectral type 


T c ft 


i It It \ 

log(L*/L Q ) 






M 




v sini 


fcrit 




(mag) 




(K) 




( M ©) 


IT) \ 

(Rq) 


(10 Mq yr ) 


(km s ) 


(km s ) 


(km s 


HD 13854 


6.49 


Bl lab 


20,260 


5.18 


(22) 


32 


0.76 


920 


97 


(328) 


HD 14134 (V520 Per) 


6.55 


B3 la 


16,300 


5.24 


24 


52 


1.45 


465 


66 


265 


HD 14818 (10 Per) 


6.26 


B2 la 


18,000 


5.40 


(30) 


52 


2.13 


565 


82 


(293) 


HD 21291 


4.25 


B9 Ia° 


10,250 


4.86 


19 


85 






30 6 


195 


HD 24398 « Per) 


2.88 


Bl lb 


20,260 


4.90 


21 


23 


0.23 


1295 


67 


396 


HD 24912 (£ Per) 


4.04 


07.5 III (n)((f)) 


35,400 


5.40 


28 


13 


0.22 


2330 


213 


550 


HD 30614 (a Cam) 


4.30 


09.5 la 


29,900 


5.86 


40 


32 


5.03 


1590 


129 


350 


HD 31327 


6.09 


B2 lb 


18,000 


4.76 


(20) 


25 






60 


(377) 


HD 37128 (e Ori) 


1.70 


B0 la 


28,600 


5.78 


38 


32 


3.15 


1910 


91 


362 


HD 37742 (£ Ori A) 


1.7 


09.7 lb 


29,250 


5.40 


25 


20 


1.05 


1860 


124 


421 


HD 38771 (« Ori) 


2.04 


B0.5 la 


23,100 


5.50 


35 


35 


1.37 


1525 


83 


379 


HD 41117 (x 2 Ori) 


4.64 


B2 la 


18,000 


5.40 


32 


52 


2.36 


510 


72 


305 


HD 42087 (3 Gem) 


5.76 


B2.5 lb 


17,150 


4.72 


(19) 


26 


0.20 


735 


71 


(359) 


HD 43384 (9 Gem) 


6.29 


B3 lab 


16,300 


4.88 


(22) 


35 


0.30 


760 


59 


(331) 


HD 47240 


6.18 


Bl lb 


20,260 


4.90 


(21) 


23 


0.31 


960 


103 


(396) 


HD 52382 


6.51 


Bl lb 


20,260 


4.90 


(21) 


23 


0.33 


900 


74 


(396) 


HD 53138 (o 2 CMa) 


3.00 


B3 la 


16,300 


5.24 


23 


52 


0.78 


865 


58 


258 


HD 58350 (r) CMa) 


2.40 


B5 la 


13,700 


5.02 


20 


58 


0.91 


320 


50 


238 


HD 91316 (p Leo) 


3.84 


Bl lab 


20,260 


5.18 


22 


32 


0.63 


1110 


75 


328 


HD 119608 


7.50 


Bl lb 


20,260 


4.90 


(21) 


23 


0.34 


880 


74 


(396) 


HD 151804 


5.24 


08 Iaf 


33,500 


6.06 


(60) 


32 


11.3 


1445 


104 


(416) 


HD 152236 (C 1 Sco) 


4.77 


B0.5 Ia+ 


19,700 c 


6.05 c 


(35) 


91 c 


5.36 


390 


74 


(99) 



a From the SIMBAD database. 

b From Abt, Lcvato & Grosso (2002). 

c From Rivinius et al. (1997). 



servations of Mason et al. (1998) exclude the existence in 
HD 24912, HD 30614 and HD 151804 of bright (Am <, 3 
mag), close companions (0.035" p ^ 1.5"). Apart from 
HD 37742, a search in the Hipparcos database reveals only 
one star in our sample with a close visual component (HD 
42087: AH P =2.58 mag and p=0.6"). 

The IRAF't tasks were used to carry out standard re- 
duction procedures (i.e., bias subtraction, flat-fielding, ex- 
traction of the spectra, and wavelength calibration) . Spec- 
tra of FeNe lamps were taken immediately before and/or 
after the stellar exposure. The spectra were subsequently 
put in the heliocentric velocity frame. The Ha line profiles 
were found to be significantly affected by telluric lines, as 
the observations were often carried out under unfavourable 
atmospheric conditions (e.g., high humidity). Therefore, we 
have obtained spectra of the rapid rotator HD 65810 (Al V) 
to create a library of telluric lines that would be used on the 
target stars. Since this procedure proved unsatisfactory, we 



t IRAF is distributed by the National Optical Astronomy Obser- 
vatories, operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the Na- 
tional Science Foundation. 



used instead the high-resolution telluric spectrum of Hinkle 
et al. (2000) degraded to our resolution. The water vapour 
lines in the pseudo-continuum regions around Ha were inter- 
actively shifted and scaled with the IRAF task telluric until 
the residuals between the observed spectra and the template 
telluric spectrum were minimized. The spectra were finally 
continuum-rectified by fitting a low-order spline3 polyno- 
mial to fixed line-free spectral regions on both sides of Ha. 
We generally obtained three consecutive exposures to allow 
a more robust rectification of the spectra and filtering of cos- 
mic ray events. Any individual exposure with evidence for 
an imperfect continuum normalization was excluded during 
the co-adding operation. The typical signal-to-noise ratio 
(S/N) at Ha in the combined spectrum is quoted for each 
star in Table 3, and lies in the range 125-335 per pixel in 
the continuum. Large velocity variations of the photospheric 
He I A6678 line, which might arise from binary motion, were 
found in HD 13854 and HD 47240 (cr~15 and 30 km s _1 , 
respectively). The spectra were hence realigned and put in 
the same reference frame. The remaining stars may have 
companions (e.g., HD 37742), but the radial velocity varia- 
tions are small, i.e., comparable to the typical uncertainties 
in the wavelength calibration (~5 km s -1 ). Hence, all the 
discussed spectral variations should be intrinsic to the star. 
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Table 2. Summary of period search (all periods in days). V mln and P maJi : lower and upper limits on the stellar rotational period 
(determined from the critical and projected rotational velocities, along with an estimate of the stellar radius); P(phot): periodicities in 
the Hipparcos data (var: stars displaying photometric variability but without evidence for periodicity; non-var: stars without significant 
variability). The light curves are shown in Fig. 2 (and in Fig. 5 for HD 14134).; V (literature): Literature values for the periodicities in 
Ha (Kapcr et al. 1997, 1998; de Jong et al. 2001; Markova 2002). The numbers in italics refer to the periodicities in the UV wind line 
profiles (Kaper et al. 1999; de Jong et al. 2001). Multiple values refer to periods derived at different epochs. Numbers in brackets indicate 
that the periodicity found is longer that the time span of the IUE observations. For the Ha periodicities derived from our observations, 
■p(Ha), '2D' and 'EW indicate that the period is present in the Ha pixel-to-pixel or EW time-series, respectively. 



Name 




V ■ 

' mm 


/ max 


I ^pnu l y 


'PHitpratiirp^ 

/ 1 11 Ll^l 111 L, 1 




HD 


13854 




(4.9) 


16.5 


5.644 




1.047 (EW) 


rilJ 


14134 


'V520 Per) 


10 


40.1 


12.823 




22.2 (2D), 12.5 (2D and EW) 


TT I~v 

rilJ 


14818 


10 Per) 


(8.9) 


31.8 


2.754 




No period search performed 


T T I~v 

rilJ 


21291 




22 


144 


26.76 


[1,2,3] 


No period search performed 


T T T~v 

rilJ 


24398 


C Per) 


2.9 


17.3 


non-var 




No periods found 


rilJ 


24912 


£ Per) 


1.2 


3.17 


non-var 


2.1±0.1, 1.96±0.11, 2.086 [4,5,6] 


2.197 ? (2D and EW) 














1.9±0.3, 1.0±0.1, 2.0±0.2 




HD 


30614 


a Cam) 


4.6 


12.5 


var 


~5.6, ~7, ~10 [5,6,7,8,9] 


No periods found 


HD 


31327 




(3.3) 


20.8 


non-var 




No periods found 


HD 


37128 


e Ori) 


4.4 


17.6 


var 


[8,10] 


18.2, 0.781 (both 2D) 


HD 


37742 


C Ori A) 


2.3 


7.97 


1.407 


~6 [5,6,8] 


^80 (2D), 50 (2D), 13.3 (EW) 














(6.3±2.1), (6.1±2.7) 


1.136 (2D), 0.877 (2D) 


HD 


38771 


k Ori) 


4.7 


21.4 


var 


[7,8,11] 


4.76, 1.047 (both 2D) 


HD 


41117 


X 2 Ori) 


8.6 


36.3 


2.869 


[1,8,12,13] 


^80 (2D), 40 (2D and EW), 0.957 (2D), 


HD 


42087 


3 Gem) 


(3.7) 


18.5 


6.807 




25 (2D and EW) 


HD 


43384 


9 Gem) 


(5.3) 


29.7 


13.7 




No periods found 


HD 


47240 




(2.9) 


11.3 


2.7424 




No periods found 


HD 


52382 




(2.9) 


15.7 


var 




No periods found 


HD 


53138 


o 2 CMa) 


10 


45.7 


var 


[8,14] 


No periods found 


HD 


58350 


■q CMa) 


12 


58.2 


var 


[8] 


No periods found 


HD 


91316 


p Leo) 


4.9 


21.3 


var c 


[12,13,15] 


No periods found 


HD 


119608 




(2.9) 


15.7 


non-var 




No period search performed 


HD 


151804 




(3.9) 


15.5 


non-var d 


7.3 [6,16] 


No period search performed 


HD 


152236 (C 1 Sco) 


(46) 


62.3 


var e 


[17.18] 


No period search performed 



a The numbers in square brackets refer to previous investigations of Ha variability in our program stars: [1] Roscndhal (1973); [2] 
Denizman & Hack (1988); [3] Zcinalov & Rzacv (1990); [4] de Jong et al. (2001); [5] Kapcr ct al. (1997); [6] Kapcr ct al. (1998); [7] 
Ebbets (1980); [8] Ebbets (1982); [9] Markova (2002); [10] Cherrington (1937); [11] Rusconi et al. (1980); [12] Undcrhill (1960); [13] 
Underhill (1961); [14] van Hclden (1972); [15] Smith & Ebbets (1981); [16] Prinja ct al. (1996); [17] Sterken & Wolf (1978); [18] 
Rivinius ct al. (1997). 

b The typical uncertainties for the quoted periods have been derived from the FWHM of the peaks in the PS and amount to: 20-30 
days (V^,50 d), 1-4 days (10<,V (d)^25), 0.01-0.2 days (1</P (d)^5), and 0.01 days (P<1 d). 
c Multi-periodic variable according to Kocn (2001). 

d Photometric variable according to Balona (1992), and references therein. 

e Irregular, multi-periodic photometric variable (Sterken, de Groot & van Genderen 1997). 



3.2 Photometry 

Coordinated B- and /-band (Johnson) photometric obser- 
vations were carried out with the 24-inch telescope of the 
Bell Observatory operated by the Western Kentucky Uni- 
versity (USA). Six stars were observed during 7 nights in 
January, February and April 2002: HD 13854, HD 14134, 
HD 14818, HD 42087, HD 47240, and HD 52382. A typical 
observation (repeated up to 3 times per night) comprises up 
to 10 short-exposure (1-10 sec) CCD images of the target 
and the surrounding field (0.53" pixel -1 scale, providing an 
4.4' x 6.6' field). 

The standard processing routines were performed with 
IRAF tasks and included bias and flat field correction, fol- 
lowed by measurements of all available comparison stars in 
the fields surrounding the targets. Pre-selection of appro- 
priate comparison stars was based on their brightness and 
non- variability. All the differential measurements performed 



on the stars from individual frames were combined into nor- 
mal observation, thus producing 1-3 data points per night 
for each star. This resulted in overall accuracy ~0.01 mag 
per data point for each comparison star. This relatively low 
accuracy was mainly dictated by the lack of appropriately 
bright comparison stars, the problem being especially acute 
for HD 14818, HD 42087, and HD 47240. For HD 42087, in 
particular, there was only one very faint comparison star in 
the field, resulting in a much lower accuracy (~0.02 mag). 
These data are not discussed in the following. 



4 OVERVIEW OF THE SURVEY AND 
QUANTITATIVE ANALYSIS 

The Ha time-series are shown in Fig.l, along with the tem- 
poral variance spectra (TVS; Fullerton et al. 1996). As can 
be seen, all stars display significant variability across most 
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Table 3. Temporal sampling of the spectroscopic observations. N: total number of spectra obtained; AT: total time span of the 
observations; Min(At), Mcan(At), and Max(At): minimum, mean, and maximum time intervals between consecutive spectra, respectively. 



Name 


Run a 


Nb. nights 


N 


<S/N> 


AT 


Min(Ai) 


Mean(At) 


Max(At) 












(d) 


(d) 


(d) 


(d) 


HD 13854 


A,B,C 


16 


19 


210 


71.1 


0.07 


3.95 


24.0 


HD 14134 (V520 Per) 


A,B,C 


17 


22 


205 


67.1 


0.07 


3.19 


24.0 


HD 14818 (10 Per) 


C 


7 


7 


170 


12.0 


0.93 


2.01 


4.06 


HD 21291 


C 


7 


7 


190 


8.01 


0.93 


1.33 


2.97 


HD 24398 « Per) 


A,B 


9 


13 


305 


38.8 


0.07 


3.24 


24.9 


HD 24912 (£ Per) 


A,B,C 


17 


23 


250 


75.8 


0.07 


3.45 


26.0 


HD 30614 (a Cam) 


A,B,C 


16 


19 


240 


70.9 


0.07 


3.94 


26.0 


HD 31327 


A,B 


10 


12 


265 


38.9 


0.10 


3.54 


23.0 


HD 37128 (e Ori) 


A,B,C,D 


17 


21 


335 


129 


0.05 


6.44 


48.9 


HD 37742 (( Ori A) 


A,B,C,D 


23 


27 


275 


133 


0.10 


5.10 


48.9 


HD 38771 (k Ori) 


A,B,C,D 


19 


22 


280 


130 


0.15 


6.17 


49.8 


HD 41117 (x 2 Ori) 


A,B,C,D 


22 


25 


220 


131 


0.11 


5.44 


49.9 


HD 42087 f3 fipml 


A B C D 


20 


20 


210 


135 
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D 


6 


6 
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a A: 2001 Nov. 18-Dcc. 04; B: 2001 Dec. 27-28; C: 2002 Jan. 21-Feb. 03; D: 2002 Mar. 24-May 07. 



of the Ha profile (the extent of the variability in velocity 
space is well correlated with the wind terminal velocity: Av 
~ ±0.3 Voo). One can note that morphological evidence for 
a Ha feature partly formed in the outflow (i.e., a line profile 
clearly departing from a pure photospheric one) is systemat- 
ically accompanied by much more prominent changes. The 
detected variability is qualitatively similar to the cases re- 
ported by Ebbets (1982). Noteworthy is the emission- like 
episode experienced by HD 14134 and HD 43384. The dra- 
matic variations previously reported in HD 21291 (Deniz- 
man & Hack 1988; Zeinalov & Rzaev 1990) and HD 91316 
(Smith & Ebbets 1981) are not seen in our data. 

Most changes take place on a daily time-scale, although 
significant hourly variations are sometimes observed. These 
subtle, short-term changes should be considered with some 
caution. In some cases, we cannot exclude the possiblity that 
they largely arise from an imperfect continuum normaliza- 
tion and/or removal of telluric features (Section 3.1). To 
better assess the reality of these short-term variations, we 
use a 'quality flag' defined as the number of co-added con- 
secutive exposures (indicated for each spectrum in Fig.l). 
Hourly changes between spectra rated as T' should be re- 
garded with some suspicion (e.g., HD 14134), whereas they 
are very likely to be real between spectra rated as '3' (e.g., 
HD 37128). Our data support the hourly changes found in 
HD 38771 by Rusconi et al. (1980), but not the variations 
on a time-scale of 30 minutes reported for HD 53138 by van 
Helden (1972). 

One should bear in mind that the observed variations in 
Ha may not be straightforwardly associated to wind activ- 
ity. It is likely that time-dependent changes, possibly arising 
from pulsations, as commonly observed in early-type super- 



giants (e.g., de Jong et al. 1999; Kaufer et al. 1997), also af- 
fect the underlying photospheric profile. Although our mod- 
est spectral resolution precludes a detailed study, we choose 
to illustrate the temporal changes in photospheric lines by 
means of He I A6678.§ With the exception of HD 152236 
which occasionally exhibits a P Cygni profile, this feature 
appears to be primarily of photospheric origin for all the 
stars in our sample. The changes affecting this line are gen- 
erally significant (see Fig.l), and are likely paralleled in the 
Ha photospheric component. Because of substantial pres- 
sure broadening in B-type supergiants, such variations will 
extend in velocity space well beyond the projected rotational 
velocity (e.g., Ryans et al. 2002). In many cases the observed 
He I A6678 variability barely exceeds the detection thresh- 
olds, which prevents us from linking the patterns observed 
in Ha and He I A6678. The limited spectral resolution also 
precludes a straightforward interpretation of the variability 
as arising either from a variable, incipient emission or from 
temporal fluctuations in the shape of the absorption profile 
(e.g., HD 31327). A direct comparison between the patterns 
of variability in He I A6678 and Ha will be presented for two 
illustrative cases in Section 6.1. 

To quantify the level of spectral and photometric vari- 
ability, we define the two activity indices a p hot and a\ pv (Ta- 
ble 4). The former was calculated both from Hipparcos data 
and from our photometric observations, and is defined as: 
a phot =cr obs -a fnstr: where CT bs and <Ti ns tr are the scatter of the 



§ We do not discuss the variations that may affect the potentially 
interesting He I A5876 line because of difficulties in defining the 
neighbouring continuum, as well as edge effects. 
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Figure 1. Ha time-series for HD 13854, HD 14134, HD 14818, HD 21291, HD 24398, HD 24912, HD 30614, and HD 31327. The spectra 
are displayed in the stellar rest frame. Consecutive spectra are shifted by 0.15, 0.10, 0.10, 0.07, 0.10, 0.10, 0.10, and 0.10 continuum 
units, respectively. The mean profile is ovcrplottcd as a dashed line. The numbers to the left-hand and right-hand sides of the upper 
panels give the number of consecutive exposures used to form the corresponding spectrum and the mean Julian date of the observations 
(ifJD-2,452,000), respectively. The two bottom portions of each panel show the TVS of Ha and He I A6678 (Fullerton et al. 1996), along 
with the threshold for a significant variability at the 99.0 per cent confidence level (dashed line). 
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Figure 1. Ha time-series for HD 37128, HD 37742, HD 38771, HD 41117, HD 42087, HD 43384, HD 47240, and HD 52382. Consecutive 
spectra are shifted by 0.12, 0.13, 0.08, 0.20, 0.10, 0.12, 0.15, and 0.10 continuum units, respectively. 
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Figure 1. Ha time-series for HD 53138, HD 58350, HD 91316, and HD 119608. Consecutive spectra are shifted by 0.10 continuum units 
in all cases. 



observations and the instrumental noise, respectively (see 
Marchenko et al. 1998). The spectral variability index, ai pv , 
is the fractional amplitude of the line-profile variations, i.e., 
the amplitude of the changes normalized by the strength of 
the feature (see equation [15] of Fullerton et al. 1996). This 
index was calculated both for Ha and for He I A6678. 



5 PERIOD SEARCH 

We performed a period search analysis by calculating the 
power spectra (hereafter PS) using the technique of Scargle 
(1982) on: (a) the photometric Hipparcos data (our observa- 
tions are not amenable to a period search) , (b) the pixel-to- 
pixel Ha time-series, (c) the equivalent widths (EWs) mea- 
sured within fixed wavelength intervals encompassing the 
whole of Ha (even in the case of P Cygni profiles). Despite 
the fact that the EW and pixel-to-pixel line-profile variations 
are not, strictly speaking, two independent quantities, a pe- 
riod search in the two data sets helps assessing the reality of 
the detected signals. Subsequent correction of the frequency 
spectrum by the CLEAN algorithm was performed in order 
to remove aliases and spurious features induced by the un- 
even spacing of the data in the time domain (Roberts, Lehar 

6 Dreher 1987). 



5.1 Photometric data set 

Periodicities in the range 1-27 days have been found for 9 
stars in our sample. Several stars display significant variabil- 
ity at the 0.03-0.07 mag level without, however, evidence 
for cyclical changes. A complex, multi-periodic behaviour 
cannot be ruled out in this case. As can be seen in Fig. 2, 
our photometric observations of HD 13854, HD 14134, HD 
14818, and HD 47240 support the periodic nature of the vari- 
ations suggested by the Hipparcos data. This demonstrates 
the remarkable long-term (years) coherency in the patterns 
of variability: folded with the appropriate periods, practi- 
cally all the new data points stay well within the scatter 
limits (mainly intrinsic, considering the typical accuracy of 
the Hipparcos photometry), without any adjustment of the 
initially chosen zero phases. 



5.2 Spectroscopic data set 

The significance of the peaks in the PS can be, in principle, 
rigourously estimated after normalization by the variance of 
the data (e.g., Scargle 1982). In the case of unevenly spaced 
time-series, however, the false alarm probability (hereafter 
FAP) is only indicative and cannot be used as a robust in- 
dicator of the presence of periodic signals in the data (see, 
e.g., Antokhin et al. 1995). While assessing the reality of 
the periods found in our spectroscopic data set, we therefore 
took into account additional criteria. Namely, we inspected 
the period-folded EW/line- profile data, and compared the 
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Table 4. Photometric and spectroscopic activity indices. 
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Figure 1. Ha time-series for HD 151804 and HD 152236. Con- 
secutive spectra are shifted by 0.20 and 0.40 continuum units, 
respectively. 



EW and 2D pixel-to-pixel PS (even though the line-profile 
variations are obviously not necessarily accompanied by sig- 
nificant EW changes). The existence of a cyclical pattern 
of variability is suggested for 9 stars out of the 17 included 
in the analysis. These cases are individually discussed be- 
low. No indication of a cyclical behaviour was found in HD 
24398, HD 30614, HD 31327, HD 43384, HD 52382, HD 
53138, HD 58350, and HD 91316. Kaper et al. (1998) and 
Markova (2002) report variations on a time-scale of about 
one week in HD 30614. A search for periodic modulations of 
the UV lines in HD 53138 proved inconclusive, in line with 
our non detection (Prinja et al. 2002). 

HD 13854: Significant power at i/~0.955 d" 1 is found in 
the EW data (Fig.3), but not in the Ha time-series. 

HD 14134: This star does not exhibit any evidence for 
binarity (Abt & Levy 1973). The photometric variability 
of HD 14134 is relatively well documented (Rufener & 
Bartholdi 1982; Waelkens et al. 1990; Krzesihski, Pigulski 
& Kolaczkowski 1999; Adelman, Yiice & Engin 2000), but 
the Hipparcos data provide the first clear evidence for a 
periodic pattern with P~12.823 d (Fig.2). Two periodic 
signals at i/i=0.045 and ^2=0.080 d" 1 (P~22.2 and 12.5 
d, respectively) are unambiguously detected in the 2D PS, 
with FAPjSO.l per cent (see Fig. 4). The latter frequency 
is also found in the EW PS (once again with FAP^O.l 
per cent) and is, within the errors, identical to the value 
found in the photometric data. The non-detection of v\ both 



a Photometric variable according to Balona (1992), and 
references therein. 

b These stars do not exhibit significant variability in He I A6678 
(see Fig.l). 



in the EW and in the photometric data sets leads us to 
take vi as the fundamental frequency. The two signals are 
likely to be simple harmonics (y2^1v\). As can be seen in 
Fig. 5, maximum Ha emission nearly coincides with maxi- 
mum light. We use the photometric ephemeris when folding 
the EWs CP=12.823 days and T =2,447,867.8), as this pe- 
riod is undistinguishable, but more accurate, than the value 
derived from the spectroscopic data. 

The clearly phase-locked nature of the variations is rem- 
iniscent of some early-type stars displaying evidence for a 
decentered, dipole magnetic field (e.g., Stahl et al. 1996). 
A single circular polarization spectrum of HD 14134 was 
therefore obtained on 29 July 2003 using the MuSiCoS spec- 
tropolarimeter at the 2-m Bernard Lyot telescope (Pic du 
Midi Observatory). The peak S/N was 250 (weather condi- 
tions were relatively poor). The mean profiles obtained using 
the Least-Squares Deconvolution procedure (LSD; Donati 
et al. 1997) show no significant circular polarization within 
spectral lines. The measured longitudinal magnetic field is 
45±145 G (G. A. Wade, private communication). 

HD 24912: do Jong et al. (1999) discussed periodic 
changes in photospheric lines that may arise from pulsa- 
tions ("P~3.45 hours). This star presents a strong case for 
a persistent, cyclical variability in the UV wind line-profiles 
with P~2 days (Kaper et al. 1997). From an analysis of an 
extensive data set of UV and optical observations, de Jong et 
al. (2001) reported a 2.086-d modulation of the line profiles, 
notably in Ha. In our data we find indication for a period 
tantalizingly close to, but apparently significantly different 
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Figure 2. Hipparcos light curves for stars displaying periodic 
variations. Zero phase is arbitrarily assigned to the start of the 
observations. The size of the symbols approximately corresponds 
to 2-cr error bars (if not shown). Our B- and /-band observations 
arc ovcrplotted as filled triangles and circles, respectively (the 
points arc vertically shifted by a constant value to match the 
Hipparcos data). In this case the size of the symbols corresponds 
to 4-<T error bars. 



from this value: P~2.20 d (with a typical uncertainty of 0.03 
d). There is also some evidence for long-term changes in the 
pattern of variability. While cyclical variations are suggested 
here in the red wing of Ha (Fig. 4), such changes predomi- 
nantly occured in the blue wing on October 1994 (de Jong et 
al. 2001). Although the periodic signal found in our data is 
present both in the Ha line-profile and EW time-series, we 
warn that its formal significance is low in the 2D PS (Fig. 4) , 
while the phase-locked nature of the EW variations is only 
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Figure 3. EW of the Ha line as a function of phase (data for HD 
14134 in Fig. 5). Zero phase is arbitrarily set at maximum emis- 
sion (or, alternatively, at minimum absorption). The uncertainties 
have been calculated following Chalabaev & Maillard (1983). 



marginal (Fig. 3). 

HD 37128: Prinja et al. (2002) discussed the UV line-profile 
variations in this star, but the limited time span of the IUE 
observations (~17 hours) precluded the detection of a cycli- 
cal behaviour on a rotational time-scale. Three signals are 
detected at the 99.0 per cent confidence level in the pixel- 
to-pixel Ha PS at: v x =0.055, z/ 2 =1.070, and i/ 3 =1.280 d" 1 
(Fig.4) . i>2 can be identified as the 1-d alias of v\ . Interest- 
ingly, the highest peaks in the EW PS correspond to the 
two remaining signals, although with low significance level, 
as confirmed by the poorly coherent patterns of variability. 
Our data therefore suggest the existence of two periods at 
P~18.2 and 0.78 d, but at a low confidence level. 

HD 37742: Line-profile variations on a time-scale of about 
6 days as suggested both in Ha and in UV wind lines (Kaper 
et al. 1998, 1999). Several significant peaks (all with FAP<1 
per cent) are found in the pixel-to-pixel PS at: !/i=5xl0~ 3 , 
1/2=0.020, 1/3=0.880, and i/ 4 =1.148 d" 1 (Fig.4). Only v A is 
found in the EW data set, but this signal is likely spuri- 
ous, as judged from the lack of coherent variations. Alter- 
natively, the EW variations appear reasonably phase-locked 
when folded with i/5=0.075 d -1 (Fig.3). This signal is also 
found in the pixel-to-pixel PS, but the confidence level is 
low. 
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HD 38771: Two signals are found in the 2D PS (see Fig.4): 
at 1/1=0.210 (FAP~5 per cent) and u 2 =0.955 d" 1 (FAP~1 
per cent), but there is no evidence for periodicity in the EW 
time-series. 

HD 41117: Several peaks are found in the 2D PS (all 
with FAP<;i per cent) at: ^i=5xl0~ 3 , ^2=0.025 (and its 
first harmonic at 0.050 d" 1 ), z/ 3 =1.045, and ^ 4 =l-085 d _1 
(Fig.4). The highest peaks in the EW PS are found at 
^5=0.025, ^6=0.065, and ^7=0.205 d _1 . Among these, 
and v-j can be ruled out from the lack of phase-locked EW 
variations. Despite the detection of 1/5 both in the 2D and in 
the EW PS, the poorly coherent pattern of variability casts 
serious doubt on the existence of this 40-d signal. 

HD 42087: This star does not exhibit significant radial 
velocity variations (Mayer et al. 1999). The highest peak in 
the 2D PS is found at ^i=0.040 d _1 , and appears both in 
the emission component and in the absorption trough of the 
Ha P Cygni profile (Fig.4). Although its formal significance 
is relatively low (FAP~5 per cent), the existence of a 25-d 
periodicity is supported by the EW data (Fig.3). 

HD 47240: There is a tentative indication that the changes 
in Si IV A1400 are repetitive over a 3- to 5-d time-scale 
(Prinja ct al. 2002). An ill-defined, long-term periodic signal 
is found in the EW data set at i/i=0.0075 d _1 , along with 
smaller peaks at ^2=0.105 and ^3=0.247 d _1 . The two latter 
frequencies are likely spurious, as judged from the poorly co- 
herent nature of the phased EW variations. Interestingly, v\ 
is also found in the pixel-to-pixel Ha PS with a confidence 
level close to 95 per cent (Fig.4). However, both the close 
correspondence with the total time span of the observations 
(see Table 3) and the sparsely sampled EW data (Fig.3) lead 
us to discard it from further consideration. Furthermore, the 
Ha morphology is typical of a fast-rotating star seen nearly 
cquator-on (Conti & Leep 1974; Petrenz & Puis 1996). It is 
therefore unlikely that spectral changes operating on such a 
long time-scale would be associated to rotation. Variations 
intrinsic to a circumstellar disc-like enhancement cannot be 
ruled out a priori, but they would manifest themselves as 
V/R variations of the line emission peaks, which are not ob- 
served (Fig.4). 



6 DISCUSSION 

6.1 Incidence of cyclical line-profile variations 

The results of our period search are summarized in Table 
2. Several of the listed periods must be regarded as tenta- 
tive, especially those detected either in the 2D or in the EW 
PS. Independent confirmation is clearly needed in this case. 
We also recall that HD 13854 and HD 47240 might be part 
of a close, massive binary system (Section 3.1). Since our 
survey was not designed to detect short-term periodicities, 
the found periods close to one day must also be regarded 
with some caution. They might arise from rotational mod- 
ulation only for the fastest rotators. Conversely, they might 
be associated with the pseudo-periodic ejection of 'puffs' of 



material. Such an outward propagation of density perturba- 
tions has been proposed to account for the recurrent changes 
(with a time-scale of ~20 days) in the Ha line of the B0. 5 
hypergiant HD 152236 (Rivinius et al. 1997). Changes on 
a time-scale of about one day might be compatible with 
the characteristic flow time through the Ha line-formation 
region in our program stars, provided that these putative, 
spatially-localized density enhancements follow a /3-type ve- 
locity law with an exponent close to unity. ^ Although such 
an interpretation cannot be completely ruled out, it should 
be noted that HD 152236 is considered as a candidate lu- 
minous blue variable and is likely not representative of the 
stars in our sample (Sterken et al. 1997). Pulsations can also 
operate on a 1-d time-scale, but we find no trivial relation 
between the spectroscopic and photometric periods (the lat- 
ter presumably related to pulsations, as argued below) . 

Summarising, we consider that the two stars with a rel- 
atively long-term periodic modulation seen both in the EW 
and in the line-profile data sets present the strongest cases 
for a cyclical behaviour: HD 14134 and HD 42087 (with 
■p~12.8 and 25.0 d, respectively). The phase-related pat- 
terns of variability are remarkably similar (see Fig. 6). In 
addition, the line-profile changes in Ha and He I A6678 also 
seem to be correlated in both stars, with excess emission in 
Ha being accompanied by a weakening of He I A6678. As can 
be seen in Fig. 5, this may also apply to weak metal lines such 
as CII A6578.05 and CII A6582.88 (note that variations in 
the continuum level may be invoked in the case of HD 14134 
in view of the large, phase-locked photometric variations). 
Since a significant amount of incipient wind emission is not 
expected for stars with such relatively low mass-loss rates 
(Table 1), this phenomenon may point to a link between 
the wind variability and photospheric disturbances. Higher 
resolution data are, however, needed to reveal an organized 
pattern of variability in photospheric lines, as observed in 
other B-type supergiants (Kaufer et al. 1997, 2002). 

The lack of cyclical variations for several stars in our 
sample can have a physical origin or may be ascribed to 
the limitations of the current data set (e.g., insufficient time 
sampling). Concerning the latter point, it is suggestive that 
the stars with a detected period in Ha are among those 
displaying the highest level of variability. The large ampli- 
tude of the variations exhibited by most stars suggests that 
an imperfect removal of telluric features or S/N issues are 
probably irrelevant. One may be more concerned by the fact 
that our observations were relatively few and spread out over 
several months. If we accept the idea that large-scale wind 
streams give rise to the variability, substantial variations in 
the global wind morphology might be expected on a time- 
scale of months/years (as in the case of the B0. 5 lb star HD 
64760; Fullerton et al. 1997). A lack of coherency in the pat- 
tern of variability on a monthly time-scale would thus not 
be surprising, especially for the fastest rotators in our sam- 
ple. It is unlikely in this case that a strong periodic signal 
would appear in the data. Another issue which may lead to 
a loss of periodic signal is the complex pattern of variability 

^ Here we do not consider the outward propagation of small- 
scale clumps of material resulting from the unstable nature of 
radiativcly-driven winds, as this would lead to a stochastic pattern 
of variability (e.g., Eversberg, Lepine & Moffat 1998). 
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Figure 4. Pixcl-to-pixcl power spectra of the H« time-series. The number of iterations and gain of the CLEAN algorithm were typically 
set to 2000 and 0.6, respectively. All the PS are variance-normalized to provide a uniform assessment of significance. The solid contours 
arc drawn for intensities of 12.9, 10.6, and 8.0 (this corresponds to significance levels of 99.9, 99.0, and 95.0 per cent, respectively). The 
dotted contours arc drawn for intensities of 7.0, 5.0, 3.0, and 1.0. In each diagram the frequency range is chosen to show more clearly 
the detected signals; no significant power peaks are found outside this frequency domain. For HD 24912, the dashed line indicates the 
2.086-d period found by de Jong et al. (2001). The bottom section of each panel shows the mean spectrum (in the stellar rest frame), 
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presented by Ha in B supergiants (see, e.g., Kaufer et al. 
1996; de Jong et al. 2001; Prinja et al. 2001), emphasizing 
the need for high spectral/time resolution. 



6.2 Evidence for rotational modulation? 

Detailed studies of the line-profile variations in the UV do- 
main have inferred a rather simple geometrical structure for 
the wind of some B supergiants, with generally two (or four) 
symmetrical, spatially-extended wind streams (e.g., Fuller- 
ton et al. 1997). Consequently, one generally expects the 
recurrence time-scale of the line-profile variations to be an 
integer fraction of the stellar rotational period, rather than 
the period itself. In HD 24912, for instance, the period found 
in Ha can be identified, within the uncertainties, with one 
half of the rotation period (de Jong et al. 2001). The fact 
that the periodicity found in the spectroscopic data set for 
HD 42087 is longer than the estimated maximum rotational 
period is therefore puzzling in this respect (Table 2). Apart 
from rotation, however, there is a clear indication for the 
existence of an additional (perhaps dominating) line broad- 
ening mechanism in B-type supergiants (plausibly tangential 
macroturbulence; Howarth 2004). The projected rotational 
velocities are thus likely to be grossly overestimated for the 
stars in our sample, with the result that the true rotational 
periods may be significantly longer than anticipated (Ryans 
et al. 2002). The modulation time-scale of Ha in HD 14134 
and HD 42087 is commensurate with the probable rotation 
periods, but rotational modulation cannot be firmly estab- 
lished in view of these uncertainties. For the former, the 
12.8-d period found is close to the estimate of the minimum 
rotational period (~10 days), which would imply that this 
star is a fast rotator viewed almost equator-on. However, be- 
cause in HD 14134 the vsini is typical of stars of the same 
spectral type (Howarth et al. 1997), this 12.8-d time-scale is 
more likely to be an integer fraction of the rotation period. 

The prominent changes affecting the P Cygni absorp- 
tion component of the Ha line in HD 14134 and HD 42087 



(Fig. 6) may diagnose variations in the amount of wind ma- 
terial seen projected against the stellar disk because of an 
azimuthal dependence of the mass-loss rate. Note in this re- 
spect that in HD 14134 maximum Ha emission and max- 
imum light do not seem to exactly coincide (see Fig. 5). 
The existence of such a phase offset (which needs to be 
confirmed thanks to a more accurate ephemeris) would be 
consistent with a picture in which photospheric perturba- 
tions lead to the development of large-scale, spiral-shaped 
wind structures (Cranmer & Owocki 1996). Let us suppose 
a bright photospheric region causing a local enhancement 
of the mass-loss rate. One observes a classical P Cygni pro- 
file when this 'overloaded' portion of the wind is projected 
against the stellar disk. However, the motion of this wind 
feature towards the stellar limb because of rotation will 
cause the growth of Ha emission. A small phase lag (i.e., 
A0<O.25) would be expected if the spot has significant az- 
imuthal extent. A substantial modulation of the local mass- 
loss rate may also be invoked to explain the large, phase- 
locked EWs variations seen in other stars (e.g., HD 13854; 
Fig.3). 



6.3 Relation to photospheric variability 

The hypothesis that the development of azimuthally struc- 
tured outflows in early-type stars is intimately linked to pro- 
cesses taking place at the photosphere (Cranmer & Owocki 
1996) has received considerable observational support (e.g., 
Henrichs, Kaper & Nichols 1994; Reid & Howarth 1996). 
Much more controversial, however, is the nature of the 'seed' 
perturbations, with pulsational and magnetic activity as the 
most straightforward candidates. In the light of our results, 
we discuss below the relevance of these two scenarios in B- 
type supergiants. 
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Figure 6. Grey-scale plots of the residuals of Ho and He I A6678, as a function of phase, for HD 14134 (P=12.823 d) and HD 42087 
(V=25.0 d). The residual profiles are the individual spectra binned to 0.05-phase resolution minus the mean spectrum. Excess emission 
(or, alternatively, weakening of absorption) appears brighter in these plots. We use the photometric ephemeris in the case of HD 14134, 
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binned profiles. The TVS, along with the horizontal dashed line indicating the 99.0 per cent confidence level for significant variability, 
are displayed in the lower portion of each panel. The grey-scale plots arc displayed in the stellar rest frame. 



6. 3. 1 Pulsations 

The Hipparcos data have revealed the existence of a large 
population of B-type supergiants with daily, quasi-periodic 
light variations. These stars might constitute the extension 
to higher luminosities of the slowly pulsating B stars, sug- 
gesting the K-mechanism excitation of multiple, non-radial 
gravity modes arising from a metal opacity bump at T ~ 
2 x 10 5 K in stellar interiors (e.g., Waelkens et al. 1998). 
Both the amplitude and the time-scale of the photometric 
variations are compatible to what is observed in our sam- 
ple. To investigate this issue in more detail, we plot in Fig. 7 
our stars on a HR diagram, along with the theoretical insta- 
bility strips for non-radial g modes calculated without mass 
loss and rotation (Pamyatnykh 1999). The stars with a clear 
cyclical behaviour (a sample evidently biased towards the 
stars with large photometric amplitudes) do not appear to 
fall in the predicted domains of pulsational instability. The 
NLTE spectral synthesis of Kudritzki et al. (1999) suggests 
systematically higher effective temperatures and luminosi- 
ties for some stars included in our study. However, such a 
temperature offset may not be sufficient to bring these stars 
in the instability zone. Furthermore, only the late B-type 



stars in our sample seem to follow the period-luminosity 
relation for g modes (see fig. 3 of Waelkens et al. 1998). 
Massive stars are also believed to experience a variety of 
strange-mode and resonance instabilities (Kiriakidis, Frickc 
& Glatzel 1993), but there is no indication that the pho- 
tometric variations are related to this phenomenon (Fig. 7). 
The ubiquity of pulsational instability in this part of the 
HR diagram (e.g., Fullerton et al. 1996; Kaufer et al. 1997) 
leads us to consider pulsations as most likely responsible 
for the photometric variations, although it is clear from the 
above discussion that the pulsation modes have yet to be 
identified in this case. The apparent dichotomy between the 
photometric properties of O and B supergiants seen in Fig.7 
(the former exhibiting significantly lower variability ampli- 
tudes) suggests that the pulsations are of different nature 
in these 2 class of objects (e.g., higher modes are excited in 
O-type stars). 

The complete absence of the photometric periods in the 
Ha time-series (HD 14134 excluded) seems at first glance 
to argue against pulsations as driver of the wind variabil- 
ity (once again under the assumption that the light varia- 
tions can be ascribed to this mechanism). As can be seen in 
Fig.7, the stars displaying the highest level of Ha variability 
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Figure 7. Upper panel: the circles show the positions of the program stars in the HR diagram (some stars have been slightly shifted along 
the ordinate axis for the sake of clarity). The size of the symbols is directly proportional to the amplitude of the photometric variations 
in the Hipparcos data, a p hot ( sce Table 4). The squares and triangles show the position of galactic OB A- type stars (van Genderen et 
al. 1989) and blue supergiants in NGC 300 (Brcsolin et al. 2004), respectively. The size of the symbols is proportional in this case to 
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periodicities, respectively. Lower panel: positions of the program stars in the HR diagram with the size of the symbols being proportional 
to the fractional amplitude of the line-profile variations in Ho, a\ pv . The evolutionary tracks (dotted lines) have been calculated for solar 
mctallicity and for the various initial masses indicated to the right-hand side of the panels. Domains of the strange mode oscillations are 
overplottcd as solid lines, while thick lines correspond to very unstable phases (from Kiriakidis et al. 1993). The pulsational instability 
strip of non-radial g-type modes is shown as (nearly vertical) solid lines. The models have been calculated using the OPAL opacities, a 
hydrogen abundance, X=0.70, and a heavy element abundance, Z=0.02 (from Pamyatnykh 1999). 



are not preferentially found close to the zones of pulsational 
instability. This hypothesis can be investigated further by 
looking for a correlation between the amplitude of the pho- 
tometric and Ha changes, as parametrized by the activity 
indices a\ pv and a p hot (the latter calculated from Hipparcos 
data; Section 4). Although there is a hint of a positive corre- 
lation (FAP~2.3 per cent), this trend disappears when the 
4 stars displaying line-profile variations that cannot be un- 



ambiguously related to wind activity, but can be of purely 
photospheric origin, are excluded (i.e., HD 21291, HD 24398, 
HD 31327, and HD 119608). In addition, there is also no cor- 
relation between a ipv (Ha) and a\ pv (He I A6678). 

The lack of a clear statistical relationship between the 
levels of activity at the photosphere and at the base of 
the outflow (as diagnosed by the photometric changes/He I 
A6678 and Ha, respectively) suggests that the same physical 
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mechanism (i.e., presumably pulsations) is not responsible 
for the photospheric and wind activity. This conclusion re- 
lies, however, on statistical grounds and should be in prin- 
ciple examined on a star-to-star basis: here we note the ex- 
ample of HD 64760, where the same periodic signal seems to 
appear in photospheric and wind lines (Kaufer et al. 2002). 
Although detailed models are still lacking, one may also ex- 
pect a complex interplay between any pulsation signal and 
wind changes (e.g., beating of some pulsation modes). This 
is likely to result in a loose (if any) relationship between 
photometric and Ha changes, for instance. 

6. 3. 2 Magnetic fields 

A growing amount of evidence suggests the existence of weak 
magnetic fields in early-type stars, although firm detections 
are still scarce and might be restricted to some atypical 
cases. Spectropolarimetric observations of the very young 
04-6 star 1 Orionis C have revealed a relatively strong 
magnetic field, 1.1±0.1 kG (Donati et al. 2002). The exis- 
tence of a weaker field in the Bl IV star (3 Cephei (with a 
polar component of 360±30 G) has also been reported re- 
cently. This field appears to be sufficient to control the out- 
flow at large distances from the star (out to ~9i?*; Donati 
et al. 2001). Magnetic activity in OB stars is also hinted 
at by X-ray observations of the 09.5 II star HD 36486 (S 
Ori) and HD 37742 (Waldron & Cassinelli 2000; Miller et 
al. 2002) which challenge the generally accepted view that 
the X-ray emission in early-type stars solely arises from out- 
wardly propagating shocks developing as the result of the in- 
herently unstable nature of radiation-driven winds (see also 
Feigelson et al. 2002). Recent theoretical studies indicate 
that a magnetic field of a moderate strength can be pro- 
duced via dynamo action at the interface between the con- 
vective core and the outer radiative envelope, and brought 
up to the surface on a time-scale typically smaller than the 
main-sequence lifetime (Charbonneau & MacGregor 2001; 
MacGregor & Cassinelli 2003). 

Apart from our observations of HD 14134 (Section 5.2), 
HD 24912 is the only star in our sample for which a mag- 
netic field has been searched for via spectropolarimetric ob- 
servations. A 3-cr upper limit of 210 G in the longitudinal 
component was derived (de Jong et al. 2001). It is worth 
noting that this does not set stringent constraints on the 
potential role played by magnetic activity, as even relatively 
weak fields (5^100 G) are suspected to significantly perturb 
line-driven flows (Owocki 1994). This issue has been recently 
worked out in more detail by ud-Doula & Owocki (2002) who 
studied the influence of a large-scale magnetic field with a 
dipole configuration on the radiatively-driven outflow of an 
early-type star. They found that the disturbance imposed 
by the magnetic forces can be parametrized by a confine- 
ment parameter, rj—B^Rl/Mvac, where B cq is the equato- 
rial field strength at the stellar surface. The magnetic field 
appears to significantly impact on the wind morphology and 
dynamics even for values as low as r)=0.1, while the emer- 
gent flow is channeled along closed magnetic loops for rf£ 1. 
For the physical parameters typical of our program stars 
(e.g., HD 14134; see Table 1), we obtain 77=1 for B cq ~70 
G. In this case the Alfven radius would be about 2 R+, i.e., 
comparable to the extent of the Ha line-formation region 
(ud-Doula & Owocki 2002). This rough estimate shows that 



magnetic activity cannot be ruled out as driver of the wind 
activity at this stage. A much stronger field is unlikely in 
HD 14134, considering our spectropolarimetric observations 
together with the lack of strong X-ray emission (Lx[0.2- 
4.0 keV] & 33.0 ergs s" 1 ; Grillo et al. 1992). Fields of the 
same magnitude can also be expected in other stars from 
our sample in view of their relatively low X-ray luminosities 
(Berghofer, Schmitt & Cassinelli 1996). 



7 CONCLUSION 

Our long-term monitoring has confirmed the ubiquity of Ha 
line-profile variations in OB supergiants. The changes result 
from the combined effect of a varying amount of incipient 
wind emission and an intrinsically variable underlying pho- 
tospheric profile, with the relative contribution of these two 
phenomena varying from one star to another. As a conse- 
quence of the large amplitude of the Ha variations observed, 
a precise determination of the stellar atmospheric parame- 
ters will not be generally possible from profile fitting of the 
Balmer lines in snapshot spectra. Taken at face value, for 
instance, the EW variations shown in Fig. 3 would translate 
in the most extreme cases into a ±0.8 dex uncertainty in 
the derived mass-loss rate (Puis et al. 1996; equation [43]). 
Although it remains to be seen whether a sophisticated mod- 
elling would result in drastically different figures, the large 
Ha changes intrinsic to the early-type supergiants may be 
one of the most important factors limiting the accuracy with 
which these objects can be used as extragalactic distance in- 
dicators (see discussion in Kudritzki 1999). 

This survey has revealed the existence of cyclical Ha 
line-profile variations in HD 14134 and HD 42087. The low 
detection rate (2 cases out of 17) may have several ex- 
planations, including observational limitations or a loss of 
coherency in the patterns of variability. HD 14134 is of 
particular interest, as it is one of the very few early-type 
stars clearly exhibiting the same pattern of variability in 
the photometric and spectroscopic modes. The existence of 
azimuthally structured outflows in our program stars might 
explain the dramatic Ha changes, but firm conclusions must 
await detailed modelling (see Harries 2000 for the first effort 
in this direction). Sensitive searches for close companions 
(e.g., via long-baseline optical interferometric observations) 
are also needed to firmly rule out a variability arising from 
wind interaction. 

We do not find evidence for a direct coupling between 
the wind activity and base perturbations arising from pul- 
sations, but we caution that our data are not well suited to 
examine this aspect in detail. It is conceivable that in sev- 
eral instances the globally anisotropic nature of the outflow 
primarily results from the putative existence of small-scale 
magnetic loops anchored at the stellar photosphere (or al- 
ternatively a presumably fossil, large-scale dipole field). A 
major observational effort is currently undertaken to detect 
such weak magnetic fields in early- type stars (e.g., Henrichs 
et al. 2003). Firm detections can be expected in the forsee- 
able future. They would place the above conclusions on a 
quantitative footing and directly address the relevance of 
this scenario. 
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